Collective excitations and low temperature transport properties of bismuth 
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We examine the influence of collective excitations on the transport properties (resistivity, magneto- 
optical conductivity) for semimetals, focussing on the case of bismuth. We show, using an RPA 
approximation, that the properties of the system are drastically affected by the presence of an 
acoustic plasmon mode, consequence of the presence of two types of carriers (electrons and holes) in 
this system. We find a crossover temperature T* separating two different regimes of transport. At 
high temperatures T > T* we show that Baber scattering explains quantitatively the DC resistivity 
experiments, while at low temperatures T < T* interactions of the carriers with this collective 
mode lead to a behavior of the resistivity. We examine other consequences of the presence of this 
mode, and in particular predict a two plasmon edge feature in the magneto-optical conductivity. 
We compare our results with the experimental findings on bismuth. We discuss the limitations and 
extensions of our results beyond the RPA approximation, and examine the case of other semimetals 
such as graphite or IT — TiS 62- 
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I. INTRODUCTION 

Bismuth is a material which plays an important role 
in solid state physics. Due to an extremely small Fermi 
surface this material provides the remarkable possibil- 
ity to observe strong effects induced by the presence of 
external fields, pressure and temperature, even if these 
external forces are of moderate amplitude. Several im- 
portant phenomena, Shubnikov-de Haas^ and de Haas- 
van Alphert^ effects, were observed for the first time in 
bismuth. In a last few years a series of experiments have 
once again drawn the attention of the community to el- 
emental bismuth, and challenged our understanding of 
this material. High pressure optical spectroscopy mea- 
surements indicate that the mechanism of semimetal to 
semiconductor decay is not fully understood^. Reflec- 
tivity measurements^ showed large changes in the plas- 
mon frequency and anomalous mid-infrared features, in- 
dicating strong scattering of the electronic degrees of 
freedom by a plasmon collective mode. An extremely 
strong Nernst signal with unusual temperature depen- 
dence (both at low^ and high fields^ii) was also reported. 

These recent experiments were showing clearly that 
the question of transport in bismuth was still not un- 
derstood. In fact similar questions still existed for the 
standard resistivity as well. The majority of the resistiv- 
ity experiments had been done in late 70s, beginning with 
the works of Hartmann^, Kukkonen- and later studie&i^, 
who showed that down to 4K the Fermi liquid theory 
(with components of very different masses) works quite 
well. The resistivity behavior at lowest temperatures 
was explained within this theory2iii. However the dis- 
cussion was not closed because one year later more de- 
tailed, lower temperature data by Uher— was published, 
and showed a significant deviation from expectations: the 
behavior changes smoothly into a behavior at the 
lowest temperatures. A quite complex theory involving 
coupling to particular group of phonons was proposed as 
an explanation of this result^. 

The purpose of this work is to re-examine the theory of 



transport in semimetals. We show that these anomalous 
transport properties have a simple explanation. They 
come from the fact that the coupling between electrons 
and holes with very different masses induces many body 
corrections to the Fermi liquid picture. In fact the above 
mentioned change in the resistivity was the first example 
how interactions in semimetals modify the simple Fermi 
liquid picture. Although our study is mostly focused on 
bismuth, we also examine other materials which have 
been recently the subject of intensive studies, such as 
graphitei4.i5' and IT - TiS'ea^JI. 

The structure of this paper is as follows. In Sec. HIl we 
introduce the model of interaction between the electrons 
and holes. We show that at low temperatures, a collec- 
tive acoustic plasmon mode exists and plays a central 
role in the properties of the material. We discuss also 
the high temperature regime of conductivity and show 
that the Baber mechanismi^ is the dominant source of 
resistivity in this regime. We then examine in Sec. IIIII 
the low temperature regime for the resistivity. We de- 
velop an effective theory for this regime and derive a new 
p(T) dependence, which is a direct consequence of the 
existence of a collective acoustic plasmon mode. We ex- 
amine the magneto-transport in Sec. IIVI We show in 
particular that a double plasma edge must exist. Some 
discussion on the validity of the approximations used to 
derive the above mentioned results are indicated in SecFVl 
and conclusions in Sec. I VII Some technical details can 
be found in the appendices. 



II. MECHANISM OF RESISTIVITY 

A. Band structure and hamiltonian 

The peculiarity of bismuth comes from the very small 
characteristic energy scales of its Fermi liquid (see Fig.[T|). 
This stems from a slightly distorted cubic crystal struc- 
ture (the distortion angle is smaller than 3°). In the 
absence of distortion bismuth would be a band insula- 
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FIG. 1. A sketch of the band structure of bismuth: three 
electron pockets and a central hole pocket, (top) the energy 
dispersion. Energies are in meV. (bottom) The four pockets. 
The scale is grossly exaggerated. In fact kp are 10^ smaller 
than intra-pocket distances. The electron pockets are slightly 
tilted out of bisectrix-binary plane Ci — C2 perpendicular to 
the trigonal axis C3. Figure after Refl20l 



tor, instead it is a quite rare rhombohedral space group 
A7 (without inversion symmetry). Bismuth becomes a 
semimetal with very small amounts of fermions active at 
accessible energies. 

The hamiltonian of those carriers close to Fermi energy 
in bismuth reads in general: 

3 

^ = ^0 + E ^O-" + H'lnt' + Htn^ + Htn^ (1) 
v=l 

In the above h denotes a hole pocket, while there are 
three electron pockets denoted by ev. The first two terms 
are the free fermion kinetic energies. We approximate the 
kinetic energy of each type of carrier by a a free dispersion 



TABLE I. The band parameters of bismuth according to 
Ref. [l^. Masses are in units of the mass of the electron, 
and wavevectors in units of reciprocal length 1.386A~^. The 
values are given along the principal axes of the ellipsoids as 
defined in Fig. [1] 



relation 

rw = ET4-£^F (2) 

where a denotes the species, the index i are the three 
principal axes of the energy ellipsoid (see Fig. [1]), mf is 
the mass tensor, the k are the momenta centered on the 
corresponding pocket and Ep the Fermi energy of the 
corresponding species. The values of the parameters are 
summarized in Table HI We want to emphasize that this 
peculiar band structure was computed quite accurately 
and this result was later confirmed by many experimental 
probes. 

The latter terms in ([1]) are interactions, whose influ- 
ence is the subject of this study. We divided interac- 
tions into three groups: hole-hole, electron-electron and 
electron-hole. Each one of these terms has the form 

HtJ = Y.V-(^{q)p^{q)pp{~q) (3) 

where a, (3 runs among the species. The density operators 
Pa{q) is given by 

^^(g) =E^i,/c+qCa,fc (4) 

k 

where the , c are the standard fermionic creation and 
destruction operators and a summation overs the spin 
degrees of freedom is implicit. The interaction potential 
is the long-range Coulomb potential 

where the hole charge is the opposite of the electronic 
one Ch = —Be and Coo is the dielectric constant due to 
the rest of the material. Because of the very small size of 
the pockets and the large distance in momentum space 
among them all interactions which implies a transfer of 
particle from one pocket to the other must occur with a 
large momentum transfer. These terms are thus poten- 
tially strongly suppressed compared to the intra-pocket 
interactions given the long-range nature of the Coulomb 
potential ([S]). 
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B. Screened Coulomb interaction: acoustic 
plasmon 



FIG. 2. The RPA series for effective interactions for the case 
of a two component system. We denote generically these 
components by r (rapid or fast) and s (slow). In the above 
diagram, the thick coil is the screened effective interaction 
14ff (5,01), the solid blue denotes using the standard diagram- 
matic convention a fast particle, while the doted red line is a 
slow particle. The black coil is the bare Coulomb potential 



of the electron-electron potential 



The smallness of the bismuth Fermi surface shown in 'OCCCOX 
FigUand incorporated in the first two terms of the hamil- 

tonian ([1]) has certain non-trivial implications for the /7ww^^ ' " " ^tww^ 

transport properties. When kp^ w lO^A the interac- 'OCGCO\,/00000' + mm<[ ^:'mm<'/j,y(nm + ... 

tion with angstrom size impurities is strongly reduced. 
The probability of intra-pocket umklapp scattering is 
also very low. Since the Debye temperature for Bi is 
around 150K it implies that at T ~ IK phonon influ- 
ence will be rather weak. Moreover thermal transport is 
proven experimentally to be ballistic. No evidence of any 
additional order parameter, generating quasi-particles 
on which carriers could scatter was found in the many 
decades of careful studies of this element. Thus in or- 
der to understand the transport properties we can re- 
strict ourselves to the Fermi liquid part of the problem 
described by the hamiltonian ([l}. 

Given the above points, the main type of scatter- 
ing entering resistivity should be the so-called Baber 
scattering^-, coming from the presence of several types 
of carrier of different masses. In contrast with all the 
above mentioned scattering processes Baber scattering is 
strongly favored by the band structure of Bi - the ratio 
of electrons and holes masses being in some directions 
more than a factor of 10. The resistivity in Bismuth is 
thus linked directly to the electron-hole interaction term 
in dH). 

Although the bare form of such a term is given by 
([3]), the presence of the two species of carriers leads to a 
strong renormalization of the bare Coulomb interaction 
due to dynamical screening. We thus need to examine the 
effects of screening on the interaction potential V°'^{q). 
Normally we have a tensor structure for the screened in- 
teraction (or alternatively the dielectric constant), but if 
we assume that the longitudinal and transverse modes 
do not mix, which is a good approximation in Bismuth, 
then we can use the standard Random Phase Approxima- 
tion (RPA) for the Coulomb potential. In the presence 
of two types of carrier the RPA equations are indicated 
in Fig. H 

The equation of Fig. [5] can be solved doing first the 
re-summation on one of the species: 



l-V--{q,Uj)U,,{q,Uj) 



(8) 



where nee(9,w) denotes the sum over the three electron 
pockets of the corresponding retarded density-density 
correlation Ilee{q,oj) = J2t=i^ee,iy{q,^)- The expres- 
sion ([U can be rewritten as 



i-Ko'ui('z)[nw.(<z,^) + nee(g,w)] 



(9) 



where in a similar way than for the holes one can ap- 
proximate the correlation function by its free (or Fermi 
liquid) value. 

As can be seen from ^ a strong resonance in V^g{q, ui) 
exists when the denominator is small, which indicated the 
presence of a collective mode. As usual the dispersion 
relation of the mode is given by 



1 = V,lliiq)Re[nhhiq,u;) + Tleeiq,Lo)] 



(10) 



y-(<z,^) 



Koul(9) 



l-V^^^Jq)U,,Hiq,uj) 



(6) 



where Uhh (q, w) is the retarded density-density corre- 
lation function for the holes, including the interactions 
which are not already present in the RPA chain of bub- 
bles of Fig. [2] The simplest and usual approximation 
for Uhh consists in taking the free correlation 11°^ (5, w). 
Using the Hamiltonian ([T]) one has 

"''"^'^'"^"fi^c. + ^(fc)-^(fc + g) + z5 (7) 

where D, is the volume of the system, / the Fermi factors, 
(5 = 0"*". Using this expression one can get the screening 



where Re denotes the real part. The imaginary part gives 
the damping of the mode. As is well known^^"^'^ for the 
case of a two component system (ITUl) has two solutions. 
One is the standard optical plasmon, but a second one is 
an acoustic mode 

^ac = Uq= (Cac - iT^^)q 

where i?e[it] = u' = Cac is the acoustic plasmon veloc- 
ity and Im[u] — u" — — r^^ the damping of the mode. 
For the case of isotropic Fermi surfaces with the same 
Fermi wavevector for the two species (neutrality condi- 
tion) and one heavy rus and light masses rrir (for slow 
and rapid carriers respectively) such that {Vps/VFrY ^ 
1 <^ ms/rrir, where the Vpi {i = r, s) are the respective 
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Fermi velocities kpi/mi with Ep = kp^/{2mi) in the ze- 
roth order approximation (extreme mass difference) , one 
has 



and 



(tO)-i = 7:Vfs/12 



(11) 



(12) 



Physically such an acoustic plasmon comes from the fact 
that for two different velocities, one of the species (the 
light one) will have a much larger Fermi velocity than the 
acoustic plasmon mode. In that case the corresponding 
Ilrr{q,i^) essentially tends to a constant HrriQ 0,i-o — 
0) which is the density of states at the Fermi level. This 
mode is fast enough to screen the Coulomb interaction. 
© would just become 



V^iiO = Cacg ^ 0) = 



-UrriuJ = Cacg ^ 0) 



(13) 



transforming the original long range interaction into a 
short range one. Using the well known result for the 
Lindhard function one gets 



Urrii^ = Cac<7 ^ 0) = 



+ 1- 



27r 



(14) 



For the heavy mode one is in the opposite limit Cac 3> Vps 
for which the Lindhard function is approximately 



=g^0) = 



TTlsKps 



'^Fs 



(15) 



Substituting into ([5]) one gets the results dTTI) and 

Such collective excitation has been already investigated 
in several systems2ii^i2i but is normally very elusive. It 
exists usually very close to the single-particle excitation 
spectrum, so it can be easily Landau overdamped. There 
were proposala^^ suggesting its measurability in artificial 
ID systems (quantum wires), but even then its intensity 
was shown to be much weaker than optical plasmon and 
was easily suppressed by disorder. However we claim in 
the present paper that such a collective mode plays an 
important role in Bismuth. Computing precisely its pa- 
rameter is a singularly complicated calculation given the 
complexity of the Fermi surface and the presence of the 
various masses. Nevertheless we give in Appendix |X] ar- 
guments for the existence of such a mode. Mostly we 
will proceed along the lines that such a mode exists and 
explore the consequences for the various transport prop- 
erties. 

Note that the mode never exists as a perfectly sharp 
mode, because of the damping ([T^ which always gives 
a finite width of the order of to this acoustic plas- 
mon resonance. In addition the mode will completely 
disappear when its dispersion relation will enter the the 
particle-hole continuum of excitations of the heavy (slow) 



species. Beyond this point the mode is severely damped 
and does not exist as a collective excitation anymore. 
The corresponding wavevector q* can be determined by 
matching the energy of the acoustic plasmon w = Cac9, 
taking into account a width of order T^^q with the edge 
of the particle-hole spectrum q{q + 2kph)/2mh. An esti- 
mate of this wavevector is thus 



q* ~ 2TO,,[Cac - Vph - T ^] 



(16) 



This condition which is essentially similar to the one for 
kc in Ref. ^23 (where Vpsq and 2kps were taken as energy 
and momentum units). 

Correspondingly q* defines an energy and temperature 
scale 



cq 



(17) 



This temperature play a similar role than the Debye 
temperature for acoustic phonons: there are no bosonic 
states to be occupied beyond. On the other hand below 
this temperature the physics of the system is affected 
by the existence of the extra collective mode ~ there are 
plasmon states whose occupation fluctuations can affect 
carriers mobilities. Based on the above we see that we 
can potentially distinguish two very different tempera- 
ture regimes: 

a. T > T* This is the high temperature (HT) 
regime. In that regime there is no collective mode and we 
can treat Bismuth as a double Fermi liquid. We can thus 
deal with the electron-hole scattering in the usual way, 
for example by solving the Boltzman transport equation, 
or equivalent approximations. We will briefiy discuss this 
rather conventional regime in Sec. Ill CI 

b. T < T* This is the low temperature (LT) regime. 
In that regime, the presence of the acoustic plasmon plays 
a central role in the interactions among particle. It will 
thus affect strongly the scattering between electrons and 
holes and lead to very different transport properties. The 
study of the consequences of such a mode is at the heart 
of the present paper and will be done for the resistivity in 
Sec. lIIII and in the subsequent sections for other transport 
coefficients. 



C. High temperature regime: Baber resistivity 

Let us now examine the resistivity itself. As discussed 
above we have to distinguish two different regimes of tem- 
perature. We examine in this section the high tempera- 
ture one, which shows relatively conventional transport 
properties, and will concentrate on the low temperature 
regime in the next section. 

At high temperature no collective mode is present and 
thus we have two well defined Fermi liquids. The resis- 
tivity is coming from the electron-hole interaction since 
this is the only term that does not conserve the total 
current. Indeed given the quadratic dispersion relation 
the current is proportional to the momentum and thus 
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the intra-species interactions conserve the total current. 
The main source of resistivity is this regime is the Baber 
scatteringii, which leads within a standard Boltzmann 
approximation to 



, m,,m.i,kf! , o T 



(18) 



where W is the e-h scattering rate and n the density of 
particles of a given species. This gives an approxima- 
tion for the coefficient A of the term in the resistivity 
p = AT^ , which is usually the quantity extracted from 
experimental data. Thus, provided that one is able to 
compute W ^ we have a parameter free fit for experiment. 
The prefactor W was computed using a Thomas-Fermi 
approximation for the Coulomb interactionM. leading to 



W=2. dq 



2k 



Fs 



VTF{qr[Hq/2kF.y 



(19) 



where kFs is the Fermi wavevector of the slower com- 
ponent. VTF('i') is the Thomas- Fermi screened Coulomb 
interaction: 



Veffiq) 



(20) 



Eoo = 100, 

rrih — 0.15 



Ref. Illl estimated the above terms using 
Klp_p = 3(4m|fci?i/7r?i2)where rrie = 0.03, 
which were known from previous independent measure- 
ments. This led to A = ^nVlcmK~^ in good agreement 
with the experimental data of that time 8nflcmK~^ (the 
first reported in ReflU), U.bn^cmK-^ (in RefH) and 
also the more recent ones^^^^ 12n^lcmK~^ for measure- 
ments along the binary axis. With the same formulas 
an even better value A = \A'nrtcmK~'^ can be found 
if we use more recent values of parameters: eoo — 88, 
77ie = 0.04, m/i = 0.14. The Bii-xSbx compound, with 
X = 0.037 (n^=o.037 ~ "-^=0/3) 'was also measurec^ and 
the experimental A = 33nflcmK~^, is in agreement with 
the expectations from (|18l) , showing that the Baber scat- 
tering is indeed the good description of the the transport 
in this regime of parameters. However low temperature 
deviation from this simplistic picture were reported al- 
ready in early works (see Ref|9| and a detail discussion 
there) 

Within the framework of the Baber scattering, one 
can go beyond the traditional Baber formula ([T8)) in two 
ways. First the large temperature behavior can be com- 
puted as welP''. At high enough temperature the 
behavior of (IT5|) crosses over, at a temperature of the 
order of Ti ~ 0.2y/TFeTFh to a linear temperature de- 
pendence. Experimental data in Bismuth shows indeed 
above T = 3K a clear dependence, becoming linear 
above T — 25K. This is in reasonable agreement with 
the estimate given by the above formula which for Bis- 
muth would be Ti ~ 35K. Second, one can refine the 
calculation of the coefficient A by taking into account 
the ellipsoidal character of the Fermi surface, rather than 



using the best fit for spherical Fermi surface approxima- 
tion (me — 0.04, rrih = 0.14). The procedure is given in 
Appendix |B] for the simplified case of two pockets only. 
The ratio of resistivities in different directions x and z is 
proportional to: 



Px 
Pz 



\/{mex + mtx) 



(21) 

Given the masses of Table U we expect that resistivity 
will be largest along the trigonal axis of the crystal. The 
value of the asymmetry w 1.4 that the above formula 
gives is slightly larger than the experimentally^ found 
one w 1.15. However one has to remember that due to 
highly anisotropic ellipsoids constituting the Fermi sur- 
face the experiment requires a very precise monocrystal 
orientation. 

However the Baber scattering does not allow to un- 
derstand the low temperature regime which in Bismuth 
corresponds to T < iK. For this regime one has to in- 
voke the existence of the acoustic plasmon. This is the 
regime that we consider now. 



III. LOW TEMPERATURE REGIME 
A. Effective hamiltonian for lowest temperatures 

As we determined in the previous section, at low 
enough temperature an acoustic plasmon mode exists. 
The very existence of this mode implies that for some 
values of the frequency and momenta the effective inter- 
action T4ff((7,w) between the carriers will have a pole. 
One can thus expect this pole to dominate the low tem- 
perature transport. 

In order to analyze the consequences of such a pole, 
we derive an effective low energy Hamiltonian for which 
we consider this collective mode as a particle. This can 
be done by replacing the term corresponding to Vcniq, uj) 
(as shown in Fig. [3]) by the propagator of a particle, rep- 
resenting the bosonic fluctuations of the plasma and en- 
suring that terms such as the self-energy terms such as 
the ones of Fig. |3] are correctly reproduced. Such a model 
is given by the coupling of electrons and holes to phonon- 
like excitations 



q 



a—e,h q 



(22) 

with such an Hamiltonian the electron self energy of 
Fig. [3] would be 



1 



iVn + i^n - £.eik + q) 

(23) 

where the i/„ and a;„ are the usual Matsubara frequencies 
and D{q, icun) is the phonon propagator 

Diq,^u^,,) = -^^ (24) 
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FIG. 3. a) The self-energy of a fermion (an electron or hole) 
experiencing Ve//(g,aj) interaction (a double coil as in Fig. [21) 
during its propagation; as explained in the text the effective 
interaction Veff{q,uj) is of the resonance form. 



The corresponding term for the original Hamiltonian is 
shown in Fig. |3] and is given by 



(25) 

One can show on this term, and the other diagrams, that 
the expansion is essentially identical provided one iden- 
tify the proper "phonon" propagator and interaction ver- 
tex Mq. This can easily be done on the spectral function 
A{q,uj) — ~^lniD{q,iu!n ui + iS). We only consider 
a; > since a; < can easily be deduced from it. The 
"phonon" spectral function is 



Aph((7,a; > 0) = S{uj - ujq) 



(26) 



The effective potential 14ff(g, i^) of ^ is more difficult to 
evaluate fully. However using the simple case of two types 
of particles (fast and slow) as discussed in the Sec. IIIBI 
one has 



Veff{q,Uj) 



StT TOr \ 1 



jL3 



^l~ci^q'/uj' + nTc,,/2 (27) 



There is a finite lifetime due to the imaginary part (|T2| , 
leading to a lorentzian spectral function. If for the mo- 
ment we ignore this finite width and assimilate the reso- 
nance to a (5-function peak we have for the spectral func- 
tion 

Aplasmo„(g,C^ > 0) = ^5(W-Cac|<7|) (28) 

One can thus directly identify the two processes with 

i^q = Cac|g| 



M, 



eO 



3^^™rCac \ ^ I |l/2 
) 



(29) 



where M^" is the zeroth order approximation for the 
strength of electron-plasmon interaction (as introduced 
in Eq. [22]) . The above identification allows us to re- 
place the problem of the screened Coulomb potential, 
by a problem of electrons and holes interacting with a 
bosonic particle. Of course this particle represents the 
quantization of the acoustic plasmon. 

Let us make some comments on the validity of the iden- 
tification (l29t . Of course the real calculation of the term 
Veff{q,uj) or Mg, specially at finite temperature would 
be more complex, and will certainly affect the quantita- 
tive aspects of the identification. However we expect the 
qualitative features of the mode identification to be ro- 
bust. In particular the frequency of the acoustic plasmon 
will of course be linear in m ~ q at small q, and in the 
same way the matrix element describing the coupling to 
the particles will go as ^ q^l'^ . Note that this specific 
dependence of the coupling constant is what makes the 
difference between the coupling to this acoustic plasma 
mode and the coupling to normal acoustic phonons. This 
will have of course consequences for the temperature de- 
pendence of the resistivity that we will explore in the 
next section. The most drastic approximation that we 
made was to ignore the finite lifetime of the mode and 
to concentrate the full spectral weight in a ^-function 
peak. This approximation is not essential, and the finite 
lifetime can of course be taken into account. It would 
simply correspond to a damping of the phonon mode. It 
simplifies however the subsequent calculations and allows 
to extract the physics in a more transparent way so we 
stick to it in the reminder of this paper. On a quan- 
titative level, both the broadening of the level u"q and 
the average energy u' q are proportional to g, as discussed 
in the previous section. For Bismuth a typical order of 
the ratio v!' jv! is u" ju' ~ 0.25 so we expect such an 
approximation to be reasonably quantitative as well. 

We can now use the standard diagrammatic analysis of 
electron-phonon interaction, keeping in mind the differ- 
ence in the matrix elements, to obtain the various phys- 
ical quantities when the acoustic plasmon is playing a 
major role. In particular the self energy (|23|) is simply 
given, after summation on the Matsubara frequencies by 



6(C.,) + /(C(fc + g)) 
^-il^n +UJq - ^{k + q) 

bi-LOq)+fm + q)) 

ivn -uJq- £,{k + q) 



(30) 



Since we are interested in the low energy dissipation we 
perform the analytical continuation ivn v + i5j^, then 
the limit — > and finally extract the imaginary part of 
Sg. This last step leads to a delta function constraints 
implementing the energy conservation (see e.g. (|40|) ). 
This equation gives us a 2D surface fio of solutions go for 
q. However the above sum on q in (|30)) should be limited 
to the values for which the acoustic plasmon exists. We 
should thus only look at values of q lower than q* which 
was defined in (|16l) . For small angles only one of the 
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masses dominates (jB5[) . Combining this with (|TT|) and 
the fact that we work in a; — ^ hmit, we deduce that in a 
Umit of smah angles (nearby direction of high symmetry) 
Cac{d) ^ cos{9) which imphes qo{0) = q* cos{9). One 
has thus always qaiO) < q* . Although potentially this 
criterion can depend on the direction of q, based on above 
considerations, we can assume here that the directions of 
q where q*{0) drops to extraordinarily small values are 
rare. In particular this means that the surface fto is well 
defined. To obtain numerical responses we will take an 
isotropic criterion for the upper cutoff. 



B. New p{T) dependence 

In the low temperature regime, the transport of the 
system, which is now dominated by the presence of the 
acoustic plasmon, can thus be described by the effective 
Hamiltonian (1221) . One has of course to keep in mind 
that because of the unusual q dependence of the coupling 
to this bosonic degree of freedom ([2^ the temperature 
dependence of the resistivity will not be the one of the 
usual electron-phonon problem. 

To compute the resistivity we use the Kubo formula 
and express the resistivity as the current-current corre- 
lation function 




Re[aa:x{i^ 0)] 

2 

TIG 

hm^^o Im[{j^{Q = 0, uj)M-Q = 0, 0))] (31) 

UJ 

We are interested in the uniform response, so we set Q = 
0, we will omit x index in the following. The current is 
simply given by 



k 



X t 



(32) 



where a is the particle species, Cq, and m^a the corre- 
sponding charge and mass in the direction of movement. 
The corresponding diagram is shown on Fig. 3^ where 
the thick line and shaded triangle indicate that both 
propagator and interaction vertex are renormalized by 
Ve//((3', cj)). The average in (1311) will be computed with 
the effective hamiltonian of the problem given by (1^^ . 
It is important to note that given the topology of the 
diagrams entering the vertex corrections there is indeed 
no double counting of diagrams if one replaces the thick 
coil of Fig. m by the bosonic excitation as done in (P^ . 

The procedure is standard^^, and follows closely the 
one for phonons. The vertex correction makes it cum- 
bersome, so we give here a derivation based on the mem- 
ory function that has the advantage to take directly the 
vertex correction into account in a simpler way. The con- 
ductivity can be expressed as 



a{uj) 



M{u) 



(33) 




FIG. 4. a) The diagrammatic form of current-current corre- 
lation function. The double line is the full Green's function 
for the electrons or holes while the shaded box denotes the 
vertex correction which is shown in part b). In b) the ladder 
series for vertex correction is shown. The irreducible part con- 
sist on single plasmon exchange (purple coil line) between the 
fermions (electrons or holes). As for normal phonons given 
the momentum dependence of the coupling constant Mq and 
propagator, it is crucial to include the vertex correction to get 
the proper temperature dependence. 



where xo is the diamagnetic term, and Mitu) is defined 



M{u) 



{F;F)^~{F;F)o 
-x{0)lu 



(34) 



where the force F is F = [j, H], and the {F; F)^ denotes 
the standard retarded correlation function at frequency 
UJ. The averages are computed with the the part of the 
Hamiltonian that commutes with the current. Given that 
M = if the current commutes with H one recovers im- 
mediately in that case from (P5)) that the system is a 
perfect conductor. In particular in the high temperature 
regime it shows immediately that the sole source of resis- 
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FIG. 5. The force-force (F-F) correlation diagram. The sohd 
blue line indicates propagator of fermion (one, for example 
slower, component of the Fermi liquid) while the purple coil 
is the plasmon propagator 



We are interested in the temperature dependence of the 
resistivity. We can thus for simplicity assume an averaged 
mass over the Fermi surface, which would affect the pref- 
actor but not the temperature dependence. This allows 
to replace — >• The term df /d£^{k) constrains, at 

low temperature, to have k on the Fermi surface. One 
thus gets 



ImM (w ^ 0) = 



6nn 



+ q)) + b{-u;,)] 



6{S,{kF +q) + Wg) - (Wg -UJq) 

{Si^ikp +q)+ UJq) + S{C{kF + q)- iUq)) 

(39) 

If the temperature is small the Fermi and Bose factors 
impose qcac ^ T, and thus the argument of the d function 
simplifies to 



tivity is the electron-hole interaction and the difference 
of masses between the two species. 

For low temperatures, in the case of the Hamiltonian 
((22)) . and using the definition of the current ((32|) one 
obtains for F (for simplicity we only kept two species e 
and h) 



F^Y. ^c.M,^^cl^^^^Ck,c. [bl, + b,] (35) 

TTla 

q,k,a 



Each species gives thus a contribution to F-F correlation 
which is shown on Fig. and for Matsubara frequencies 
equals to: 



1 



132^2 ^2 



" ''G{vi,k)G{v2,k + q) 



D{ujn + vi-V2-,q) (36) 



The frequency summation can be performed which 
gives 



1 Apql 
{F; F)^ = —- y — 



[fm + q)) + b{~UJ,)] [/(C(fc) - fm + q)+ UJ,)] 
iujn + S,ik) - ^(fc + q) -UJq 

- {Ojq -UJq) (37) 

After the analytic continuation one gets for the imaginary 
part of the function M{oj) 



k.q 



mk) - ak + q)- ^q) - K ^ -w?) (38) 



5{i{kF + q) + ujq) = — gcos(0)±cac(% (40) 
m 



leading to 
ImAf (w 0) 



1 



M^q^ 



dq- 



[f{LOq)+b{LUq)] (41) 



where we applied the same reasoning as for the self- 
energy (pO)) . The temperature can be rescaled out of 
the above integral leading to 



ImM(cj ^ 0) = 



0q' ][f2q3 



dq 



-[f{uq) + b{uq)] (42) 



where the / and b and the Fermi and Bose factors with 
/3 = 1 and we took isotropic case. At low temperature 
the integral tends to a constant, which can be evaluated 
in the isotropic case, so the resistivity has a 



p(T) cx 



(43) 



temperature dependence. 

We thus see that the same electronic mechanism which 
at higher temperature was giving the conventional Baber 
behavior will smoothly lead to a behavior when 
acoustic plasmons begin to govern screening. The new 
behavior comes essentially from two approximations: 
the linear dispersion of bosons and Mq ~ The upper 
limit of the integrals in (|4T|) - is the largest possible value 
of plasmon momentum. This critical wavevector q* , at 
which Landau damping suppress the plasmon as a well 
defined particle, plays a similar role than the Brillouin 
zone boundary for acoustic phonons. As we discussed in 
the first section, one can also define a corresponding tem- 
perature T* (see (fT7|) . which is the analog of the Debye 
temperature and plays a similar role in the resistivity. 
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r* T/6 r LogT 



FIG. 6. The temperature dependence of resistivity in bis- 
muth. Three regimes are distinguished ~ T, ~ ~ T^. 
The lowest temperature regime is due to the presence of the 
collective acoustic plamon mode. The two high temperature 
regimes are due to the conventional Baber scattering mecha- 
nism. 



If we take a zeroth order approximation for the inter- 
action M° (see ([25])) and T* sa IK, we can estimate the 
A coefficient in the Deybe law (as defined in Eq. 1 of 
Ref ■ [T2h . Along the trigonal axis we get A « 0.25 [/^ilcm] 
(in general A g (0.25, 0.6) [/xficm]) which is not far from 
the experimental value A « 0. 15 [/iilcm]. The zeroth or- 
der approximation is overestimated because c^^ (and also 
c^^^)\s known to be overestimated, plus as we explained 
while evaluating the integral (j4l]) not all the 2D surface 
of solution might be present. 

A full temperature dependence of the resistivity is 
shown on Fig. [51 Note that the measured T* will be 
always slightly smaller than the RPA one T^p^^ due 
to many plasmon processes (which are not accounted 
in RPA approximation) enhancing the Landau damping 
around T*. In fact one can expect that the RPA pre- 
diction Tpp^ works better for polycrystals where acous- 
tic plasmons have more decay channels (for example 
they can scatter on surface waves) alternative to Landau 
damping. This was indeed observed experimentally^^. 

Let us finally make a comment on the anisotropy of 
the resistivity, to argue that the low temperature con- 
ductivity is not strongly affected by a particular choice 
of q direction. In order to estimate the effect we propose 
a zeroth order approximation in which ~ (-^g 
where an average is taken at a given low tempera- 
ture. The direction dependence comes from the mass 
anisotropy and in this respect we know that two factors 
behave in different ways. 

• from the reasoning given in Sec. IIII Al in particular 

while at the lowest order ^ Vps ^ {iT^sq)~^^^ , 



which implies: 

(M°),-^(m,,-)-3/4 

We expect this dependence to be stronger when we 
reach q*. 

• the average density of plasmons at the lowest 

temperatures will be determined by the Taylor ex- 
pansion of the Bose-Einstein distribution b{cac{(f)q) 
which means {b^^q) ^ rug^ while at the higher tem- 
peratures (around T k,T*) the density of available 
slow fermions (necessary to build up a collective 
excitation) sets the limit which is also ~ TOs^. 

Both contributions to scale in an opposite way with 
the mass thus we expect the overall direction dependence 
to be rather weak (especially because a more refined ap- 
proximation for p^ would need to contain some average 
over the whole Fermi surface). 

We expect a stronger influence of anisotropy on q* , be- 
cause this quantity depends only on the plasmon velocity 
u = Cac — iT~^. The problem is rather complex and we 
present a discussion of it in the Appendix [A) Consider- 
ing the anisotropy of q* bring us to an approximation we 
have made while deriving (|42p : we assumed that the plas- 
mons propagate equally well in all directions. As shown 
in Appendix |X] when more refined approximations for the 
polarizability n(fc, w) are used one gets more isotropic q* 
values. Indeed the measured anisotropy of T* is quite 
weak—. This makes the derivation of (|42p self-consistent. 

C. Electron-plasmon coupling non-perturbative 
theory 

As we saw in the previous sections, the presence of a 
collective mode has e drastic influence on transport prop- 
erties. The mapping to the effective Hamiltonian ((22|) 
allows to describe these effects. In the previous section 
we looked at the consequences of such a coupling to the 
collective bosonic mode, in a lowest order calculation in 
the self energy. 

However it is well known, for the case of electron- 
phonon coupling, that the combination of the electrons 
and the phonon degrees of freedom can be carried out 
beyond this lowest order and lead to a new composite 
quasi-particle, the polaron. For the case of an acoustic 
plasmon, a similar phenomenon occurs^^. In our case 
the corresponding quasi-particle would be a plasmaron 
(electron accompanied by plasmons) . The creation oper- 
ator of a plasmaron can be written using a variant of the 
Lang-Firsov unitary transformation: 

= 4 exp(-^) (44) 

where: 



10 



where ag ~ Mq/Ep measures the strength of the boson- 
fermion couphng (in (|45p we assumed it is real). The 
shift of fermion energy Re^hoson]) and mass 

boson 

]) are known: 
|A£:| = (a + 0.0123a^ + 0.00064a^ + (46) 

l^m{q) = {a/Q + OmAa^ + ...)m{q) (47) 

where m(g) indicates the mass in the q direction. The 
above formulas were obtained in the so-called large plas- 
maron picture. In the case of Bismuth the plasma is not 
rigid enough (contrarily to the ionic lattice) to sustain 
self-localization, so a small plasmaron picture is unphys- 
ical. 

Naturally, when a — > all the above effects disappears, 
but the question about the strength of full Mq (not M°) 
is open and quite difficult to access theoretically. We can 
only give an upper and lower limit for a. First, Mq must 
be smaller than the bare Coulomb interaction in order to 
make the plasmon a well defined particle. The ratio of 
bare interaction to kinetic energy is given by and as 
discussed in detail in Sec. IV Ai r, < 1 which implies a < 1. 
The lower bound of a can be obtained experimentally: re- 
cent optical spectroscopy measurement^ showed the exis- 
tence of an (optical) plasmaron in Bismuth which shows 
that the a coupling is sufficiently strong {a ^ 10~^) to 
give observable effects. 

The notion of plasmaron allows for some developments 
of the theory beyond RPA. The numerical evaluation is 
left for future investigations. We will come discuss some 
of these effects at the end of Sec. IV Al However a brief 
inspection of formulas (I44II47I) has several consequences: 

• The dominance of the linear term in (I46I47P con- 
firms the accuracy of the first order approximation 
for electron-boson even for a ~ 1, which validates 
our perturbative approach in Mq. 

• A plasmaron remains a well defined fermionic par- 
ticle, so the low energy model remains valid 
but with renormalized parameters. We can thus 
expect the behavior of the resistivity to be still 
obeyed even in an intermediate coupling regime. 

• from (|T7)l we see that in the presence of a plas- 
mon cloud the mass of the heavier fermions fur- 
ther increases by Ams, which should give rise to 
a q* increase. More importantly the ratio of the 
geometric series which determines Ato^ is propor- 
tional to mr/mr which means that the effect is 
the strongest in the direction where q* was smaller. 
Higher order effects thus tend to stabilize the acous- 
tic plasmon 

• by definition (j45p a fermion is mostly accompa- 
nied by plasmons propagating in the same direc- 
tion {qferm \\ qboson)', thcsc plasmons are also me- 
diating the effective interaction between electron 



and holes and thus the part of interaction \\ qferm 
is enhanced. This goes towards the unidirectional 
scenario advocated in Appendix [X] 

IV. OPTICAL SPECTROSCOPY: 
MAGNETO-PLASMON 

An obvious question is whether this new quasi-particle, 
is responsible for the lowest temperature resistivity data, 
can be probed by other measurements. 

Let us consider the consequences of the existence of 
the plasmon mode when there is a magnetic field acting 
on the system. We assume that the plasmon propaga- 
tion direction is perpendicular to the external magnetic 
field, the geometry which is usually used in optical spec- 
troscopy. We also assume that in this plane electrons 
are heavier than holes (this can change depending the 
bismuth sample orientation). As we will the low energy 
collective excitation can be probed by measuring the op- 
tical conductivity of the system. 

A. Dispersion relation 

The inconvenience of studying the acoustic plasmons 
with optical methods comes from the fact that in the 
limit q 0, which is the resonance condition with pho- 
tons, by definition we have w — >■ 0. The way to overcome 
this difficulty is to introduce a magnetic field. In the 
case of a standard (optical) plasmon the frequency of 
collective excitation (propagating perpendicular to the 
magnetic field) is increasing by the magnetic field as 

uj = '\Jojp + Wc • This phenomenon can be explained in- 
tuitively: the movement of a charge under a magnetic 
field requires an extra kinetic energy (to overcome the 
magnetic field vector potential. This result was found in 
several different ways: using RPA polarizability, equa- 
tion of motion technique or from simple hydrodynamic 
picture of fluctuations. The reason why such a remark- 
ably simple relation holds comes from Kohn's theorem: 
in a translationally invariant system, for q — >■ 0, the mag- 
netic field does not affect the inter-electron interactions 
encoded in ujp. This quite general theorem holds also for 
multi-component plasmas. 

Indeed for the two component plasma problem under 
magnetic field the following, a relation, valid for any value 
of the magnetic field, in the g — limit was found^: 

^{-Lof + ujf - ujf + uj'i^Y + Aujfuf] (48) 

When |a;^ — is small we find the intuitive relations 
uj\ « ujf+Lof+uj'^'^+uj'^'^ andwi « \ujloj^\+0[q), which 
allows for a clear identification of optical and acoustic 
plasmon when i? — >■ 0. We see that the acoustic mode 
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develops a gap, which should be observable in optics as a 
second, lower frequency plasmon edge. In the limit when 
\u)l — I ^ l^^pli I'^pl the uj± modes in (j48|) become hole 
and electron like respectively. It is relatively simple to 
reach this limit (when the second term under square root 
is small) in bismuth. The reason is that we always have 
rrir <^ ms, which implies that |a;^ — ~ 
We thus expect w_ to have an electron-like (in general 
heavier particle-like) magnetic field dependence. 

From (|48l) we can deduce how the frequency of the 
lower plasmon edge evolves with magnetic field. The 
temperature dependence of u!± is also frequently stud- 
ied in optical experiments. The temperature dependence 
enters (|48|) via ujp''^ ^ n^ji. The density of carriers ne,/i 
can change with temperature due to thermal excitation 
from the valence band. Because of the very small gap at 
the bottom of electron pockets A « 13.7meV (see Fig.[T]) 
this effect is particularly important for bismuth. Then a 
second type of holes, as light as electrons, will emerge, 
and due to the global compensation of electrons and holes 
Ue = Uh + nih the total number of carriers will increase. 
In the low field limit only the standard (high energy) 
optical plasmon frequency oj^{T) will be affected. 

In the high field limit the uj- will acquire a 
rather strong temperature dependence since uj'^{T) ss 
{u!p{T))'^ + (w^)^ (assuming that the electrons are heav- 
ier). From the charge neutrality condition the approxi- 
mate formula for ujp{T) can be giveit^: 



4e^ 



. 3/2 



-Lii/2{eMPi-EFe - ^Ji)))(kBTf'^ 

(49) 

where e is the electron charge, me is geometrical average 
of electron masses and m% is a mass in a given direction 
(all masses in electron mass unit). ii]^/2(6xp(/3(— i?i?e — 
/i))) is the value of the incomplete polylogarithm function 
of order 1/2: 



(50) 



I-ii/2{a) = / dx rxTT 

7o aexp(a;)-|-l 



This gives the temperature dependence of the w- for a 
very large magnetic field when this mode merges with 
the electron mode (in general the heavier carrier mode). 



B. Longitudinal f-sum rule 

The third quantity accessible from optical measure- 
ments is the relative intensity of two plasmon edges. This 
can be studied quantitatively using the longitudinal f- 
sum rule. 

In a one component plasma it can be proven that the 
optical plasmon is the only collective excitation as o — > 0. 
It comes from an exact sum rule (see p. 358 in Ref|28l): 



j (hjujlvi{e{q,ujy^) = --ujI 



An analog sum rule was derived for a multi-component 
plasma (see for example Ref. |32| ). It was also shown^ 
that in the presence of a non-local (q-dependent) poten- 
tial the right hand side of ([FT]) can change. This is par- 
ticularly important for a two component plasma when 
the non-local part might differentiate between intra- 
and inter-component interactions (defined as {q) in 
(HI)- ). Then the sum rule is not completely exhausted 
by the cj+ mode and there is some spectral weight left for 
the ijj- plasmon'^ . These findings prove that the finite 
spectral weight of acoustic plasmon does not violate the 
longitudinal f-sum rule. 

One can quantify this last statement by computing the 
strength of each plasmon, which in fact corresponds to 
/m[Ee((7)] evaluated in Sec. IIII Al 



W±{q) = 



\dRe[e{q,uj)]ldLo\, 



(52) 



where e{q,uj) is a dielectric function defined as the de- 
nominator of Eq. M e{q,uj) ^ 1 - V^^^^^{q)[nhh{q,uj) + 
Ilee{q,uj). Using the formulas given in Sec. |ll] and the 
Lindhard approximation for polarizability, we find (in the 
q^ limit): 



dRe[e{q 0,uj)]/duj\^=^ 



(r„ - R) ln( — ) 



- Vpsq 
(53) 

where and R are characteristic, constant parameters 
of the two component plasma (they are of order 1, for 
precise definitions see Appendix]^. Taking into account 
lims^o"^- = Cqc? (neglecting lm\u}\) we deduce that the 
numerator of the logarithm limits W- (q) : 



lim W^(q) 



[{r,~R) H{Cac-VFs)q)Y 



(54) 



(51) 



From (j54l) one recovers the Landau damping formula: 
when Cac — >■ Vps, W- — 0. The larger Cac is with respect 
to Vfs, the larger spectral weight the acoustic plasmon 
will contain. From (IMl) we also see that, when Cac Vps, 
the W-{q) is an increasing function of q; it is difficult to 
expand this formula in q — > limit, but using I'Hospital 
rules (for the derivative of ln(g)) one can confirm our 
argument given in Sec lIII Al that Vq'^ ~ q. 

The other implication of (|52)) is that inducing any gap 
in uj-{q = 0) spectrum will increase W-{q = 0). This is 
precisely what magnetic field does as shown in ([^S]). In 
particular in a strong magnetic field w_ — >■ uj^+uip which 
means that the acoustic mode becomes a mode of heav- 
ier (electron in our example) plasma component. Then, 
there must be also a redistribution of the relative spec- 
tral weights of the two plasmon edges uj± ; the a;_ mode 
will be more and more pronounced. To be more precise 
in the equation (1521) . to get the intensity of w_ mode, 
we substitute e{q,uj) of the heavier component plasma. 
The Re[e{q — >■ 0,w_(g))] — > q^ so the denominator of 
((52| gradually approaches zero and does it in the same 
way than Re[e{q — > 0,a;+(g))]. This resembles the fact 
that the spectral weight W-{q — >■ 0) must become as sig- 
nificant as W+{q — ?► 0). As the magnetic field increases 
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property low B large B 

cj_(r) T independent ~ t^p(r) 

W"- (B) = l/e(6J-) I W+ > W- ~ ^ I W^^W+ 

FIG. 7. Summary of the expected field dependence of the 
two plasmon edges. We present the magnetic and tempera- 
ture dependencies of the lower energy plasmon edge and the 
magnetic field dependence of its amplitude. The low and high 
field regimes are distinguished 

{q 0) gradually acquires all the density fluctuations 
of the heavier carriers. 



C. Summary of the expected effects 

The findings of two previous sections are summarized 
in Tab. [7l In order to define the border between low B 
and large B we define the ultraquantum (UQ) limit for 
which the magnetic field freezes the motion of carriers 
in a plane perpendicular to it. Then lOc dominates over 
Wp. Due to the particular band-structure of bismuth and 
the field dependence of the chemical potential electrons 
and holes enter to UQ regime at nearly the same field (~ 
lOT). Thus the above mentioned condition joj^ — wj | <C 
ujp is valid up to quite high fields. The limit of validity of 
this condition is for fields smaller but of the order of the 
ones of the UQ regime. The value will change depending 
whether one works with pure Bi or Bii-^Shx (xjO.O?) 
where the UQ limit is reduced. 

There are effects not included in the above analysis. 
First we observe that, because of charge neutrality re- 
quirement, and the small value of the gap below the elec- 
tron pocket, the chemical potential does depend on the 
magnetic field. Thus we expect: 

^ n{B) ^ i^ojl^^iB) (55) 

but this extra dependence is significant only for fields 
above the UQ regime. In the same regime one may also 
expect an anomalous Zeeman splitting and large magne- 
tostriction which in principle can change the band pa- 
rameters (masses). A large spin-orbit coupling in Bis- 
muth is responsible for these effects. They are beyond 
the scope of this paper, which is dedicated rather to low 
field effects, and we leave them for future investigations. 

The other source of finite gap in the spectrum of acous- 
tic plasmon at g = is the finite tunneling probability 
Tj^ between two types of carriers (electrons and holes in 
our case) . According to Ref. [1^ the additional gap IS.sas 
(using the notations of this paper) will be proportional 
to tunneling (or scattering) probability (see Eq. 26 in 

Ref. M) ^SAS - V^(9-0)Tx. 

Due to Bi band structure the energy-momentum con- 
servation strictly forbids such a scattering, thus at T ~ 
=> = ^sAS = 0. The situation could be 
slightly more complicated at finite temperature for which 



we can find a combination of energy conserving states, 
but the conservation of k is still never fulfilled (electron 
and hole pockets are very far in momentum space). We 
conclude that this type recombination processes, even in 
the higher order scattering events, should not affect the 
zeroth order prediction for uj^{B 0,T) ^ given be- 
fore. 



V. DISCUSSION 

In this section we want to discuss the limits of validity 
of our theory and the physics that can be expected when 
we reach those limit. The experimental relevance of these 
effects will be presented. 

A. Validity of RPA approximation 

In order to assert the validity of analysis of the previous 
sections, which is mostly based on an RPA approxima- 
tion of the interaction terms, it is important to estimate 
first the strength of the interactions in the case of Bis- 
muth. At ambient pressure the ratio of the potential 
energy to kinetic energy are respectively « 0.2 and 
« 1.5 for the electrons and the holes. These rather low 
values are due mostly to the very high, high frequency, 
dielectric constant of the background Soo — 88. It is also 
important to note that a part of the Coulomb interac- 
tion, namely the Hatree-Fock terms with exchange in- 
teraction Vcoui('z)/eTF('J') with Thomas-Fermi screening 
was already taken into account during a self-consistent 
band structure calculations, leading to renormalized dis- 
persions ee,h{k)^- 

Although this question is difficult to address quantita- 
tively it is useful to estimate if for the particular case of 
Bismuth we could trust the RPA approximation. Usu- 
ally the following arguments are used to justify the RPA 
re-summation shown on Fig. [3] are: i) the small num- 
ber of carriers implies that the long range character of 
interactions plays a major role {V{q) large for q 0). 
Then the diagrams with the largest divergence number 
(RPA series) are the most important. This is the case in 
bismuth despite the simultaneous validity of the "dense 
plasma" regime Ts < 1; ii) we are primarily interested in 
the limit of long wavelength excitations (below q*) thus 
the influence of local field effects is moderate. 

In addition to these two arguments that support the 
RPA approximation, the particular combination of band- 
structure parameters of bismuth contributes to a further 
suppression of the diagrams outside the RPA series. In- 
deed in Bismuth the Fermi surface consists of distant 
pockets for electrons and for holes. A significant part of 
the exchange diagrams (the inter-pocket ones) requires 
quite large momentum exchange, so in the limit of small 
momenta, as discussed above, (Vcsiq oo) — > 0) their 
contribution must be very small. In addition such pro- 
cesses are irrelevant in the RG sensed. Another con- 
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tributing factors comes from the fact that electrons along 
the trigonal axis can be considered as Dirac particles, 
while holes are very light in the perpendicular plane. 
If one take a (one component) gas of spinless fermions 
with linear dispersion then bubble approximation (each 
bubble with only two external interaction lines) is exact. 
In particular, in one dimension, for such a Dirac spec- 
trum, the RPA approximation would be exact. In higher 
dimensions one can thus expect that the amplitude of 
non-bubble diagrams should be at least be reduced when 
m — > 0. Finally both the large mass anisotropy (on the 
surface of each pocket) and the strong spin-orbit cou- 
pling (lowering the orbital momentum je quantum num- 
ber (je — -\- s"e) contribute to the weakening of the 
electron-electron exchange processes. Although none of 
this arguments is of course rigorous or final, they sug- 
gests that the RPA approximation in Bismuth is indeed 
a very good starting point to understand the properties 
of this material. 

Improving above the RPA is of course a very diffi- 
cult proposal. Among the diagrams one would have to 
consider there are multi acoustic plasmon processes (op- 
tical plasmon has a high energy so we can safely ex- 
clude it). Some of them are clearly suppressed. First, 
Mq — which implies that tadpole-like diagrams medi- 
ated by acoustic plasmons gives no contribution. Second, 
a multi-plasmon scattering processes has a natural cut- 
off q* above which they are strongly Landau damped. 
This mechanism of q*g reduction was invoked in the con- 
text of the Fig. m In fact given the acoustic nature of 
plasmon excitation one might be tempted to introduce a 
variant of the Migdal theorem in order to exclude pro- 
cesses with crossed plasmon lines. The Migdal theorem in 
our case has the following form: if we compare the avail- 
able momentum phase space we find that n-crossed multi- 
plasmon lines will generate corrections proportional to 
X'^ajq/Ep, where A ~ Mq. In order to estimate the value 
of the interaction we note that the electron plasmon cou- 
pling cannot be larger than the bare Coulomb interac- 
tions A < e^kp ^ Ep. Similarly the energy transferred 
by the acoustic plasmons is of the order of <T*< Ep. 
Given the hierarchy of energies between T* and Ep we 
see that even without the large mass ratio that exists in 
the case of the electron-phonon coupling, here we have a 
justification of a Migdal-like theorem to drop the higher 
order crossed plasmon diagrams. 

Of course, even after excluding all the above mentioned 
processes we are still left with many diagrams, which can 
be constructed around a dominating RPA series. The 
first class of multi plasmon diagrams which are left are 
the rainbow ones. Usually diagrams of this type are re- 
sponsible for renormalization of particles bare dispersion, 
which we have accounted during the discussion of the 
plasmaron in Sec. lIII Cl The second class of diagrams are 
the vertex corrections inside a single electronic bubble. 
The third class which would make the problem highly 
non trivial corresponds to interactions between several 
electronic bubbles which would lead to an interaction be- 



tween the plasmon modes or interactions among several 
bubbles of the ladder series. How to take into account all 
these terms is of course going well beyond the scope of 
the present study. Based on various arguments one can 
expect on a phenomenological level a renormalization of 
the plasmon velocity. A phenomenological model that 
could account for those extra contribution would be to 
add an interaction term between the plasmons in our low 
energy Hamiltonian ((22)) of the form 

^fpi-int ^ XI 9{q)bi^+qbi,^qh^bk2 (56) 

ki ,k2 ,q 

the interaction between plasmon term g{q) can in prin- 
ciple be estimated from the above mentioned diagrams, 
and this is left for a a future study. This phenomenologi- 
cal term could thus open the route to tackle the situation 
rs > 1. 



B. Pressure induced semimetal decay 

One way to control the physics of the acoustic plas- 
mon could be to consider the effect of pressure on these 
materials. Indeed the tiny Fermi surface in Bismuth is 
quite fragile versus Sb doping or pressure - the pock- 
ets empties and the material becomes a semiconductor. 
These dependencies are well established both on theo- 
retical (DFT pseudo-potentials method in Ref. HI) and 
experimenta l^^'^^ side. Obviously the collective excita- 
tions such as the plasmon will be affected by such a 
semimetal-semiconductor (SM-SC) transition. 

At the first naive level the change of carrier density 
affects only optical plasmon. The velocity Cac (see Ap- 
pendix |X| depends on the ratio of masses and Fermi ve- 
locities of the two types of carriers. These crucial param- 
eters nie.h or Vp seem to be constant up to very low car- 
rier concentrations'^^. However when n — ^ 0, the screen- 
ing of the interaction is drastically affected and Ts — oo. 
The corrections beyond RPA, discussed in the previous 
section, start to be important. On the other hand we 
know that, at least for Sb doping, the gap (at L points) 
is being closed, thus Soo increases, extending the validity 
of the "high density plasma" regime. It means that it 
should be still possible to define plasmons as quasiparti- 
cles of the system and the physics which we described in 
the previous sections should still be applicable. 

Non-trivial effects may arise only near the critical dop- 
ing Sc or pressure Pc for which rg — > oo. Two elements 
will play a role: i) as described in the previous section 
(and at the end of Sec. IIII Ap the most important ef- 
fect emerging with an increasing is the appearance 
of strong plasmon-plasmon interactions. The interaction 
term (|56p could potentially lead to an instability with a 
new minimum of the energy at a finite q. In that case 
the corresponding energy gain would favor a semimetal 
(correlated liquid) versus a semiconductor; ii) contrarily 
to ujaciq) = Cacq, the frequency of the optical plasmon 
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ujp ^ n is decreasing when n — >■ 0. At some point the 
two collective excitations will merge. This should change 
both their dispersion as well as the physics of the system. 

On the experimental side, the regime we are dis- 
cussing above was recently studied by means of optical 
spectroscopy^. A transfer of the spectral weight to the 
plasmaron peak was observed. This implies that, as we 
expect, plasmons are able to survive (or can be even en- 
hanced) . In these experiments a deviation from the Fermi 
liquid theory was found, with an abnormal rigidity of the 
metallic phase near pc- What was observed can be inter- 
preted as an abnormal increase of the collective mode 
mode frequency (deviation from a single mode RPA at 
the lowest temperatures). The extension of the theory 
proposed in the previous section, ([56| can perhaps be 
used in such a regime to explain these effects. 



C. Comparison with other multi-valley semimetals 

As already mentioned in the introduction, a few other 
examples of the multi- valley semimetals are known, with 
a band structure similar to Bismuth. These other sys- 
tems are thus potentially described by a similar theory 
than the one that we have introduced in this paper. 

The most intensively studied material in this category 
is graphite. In this case electrons and hole pockets are 
placed alongside, and are quite similar in shape, thus 
the carriers mass (and velocities) ratio is close to ~ 1. 
In that case (see Appendix q* shall be very small 
and the acoustic plasmons are always overdamped. We 
thus expect a standard Baber resistivity down to the 
lowest temperatures. This is indeed what is observed 
experimentally^Si. 

The second material, recently investigated, is IT — 
TiSe2- The band structure resembles strongly the one 
of Bismuth. It consist of three electron and one hole 
pocket with large mass differences^^, mh = 0.23, m^x = 
5.5,mey = 2.2 and Soo — 44 (due to the large Se 
polarizability)^. This suggests that we can have an- 
other family of materials where acoustic plasmons play 
a major role. The band structure gives a finite amount 
of acoustic plasmons in the system (k* ^ 0) but also a 
quite large Ts, which resembles more the situation de- 
scribed in Sec. IV B1 than the one of pure bismuth. This is 
most likely the reason why experimentally the low energy 
physics of IT — TiSei is different than in Bi. A density 
wave phase transition was measured at T^vi/ — and 
the related reconstruction of Fermi surface was revealed 
by ARPES measurements^®. Interestingly the unusual 
increase of resistivity'*'^ (attributed to an increase of the 
electron-hole scattering rate at the transition) just above 
the Tdvk seems to be similar to what was observed in Bi 
very close to the critical pressure Pc- The formation of an 
excitonic liquid was suggested to explain the transition 
at Tdw ■ Whether one can make contact between such an 
excitonic liquid description and the interacting plasmon 
theory that we described in the previous section is a chal- 



lenging theoretical question that we leave for future stud- 
ies. Note that at the lowest temperature in IT — TiSe-i 
superconductivity was recently observed^^ suggesting ad- 
ditional lines of studies that would use the similarity be- 
tween the plasmons and the acoustic phonons as an ex- 
change particle. 

VI. CONCLUSION 

We have presented in this paper a theory of trans- 
port in semimetals, concentrating specially on the case 
of Bismuth. We have shown that the physical properties 
of these systems are dominated by the presence of an 
acoustic plasmon mode at low temperatures. This mode 
that we derived in an RPA approximation of the interac- 
tions lead to a drastic change of the transport properties, 
in particular their temperature dependence, compared to 
the standard Baber mechanism which is normally invoked 
for such materials. We showed in particular that it would 
lead to a T^ behavior of the resistivity below a certain 
energy scale T* dependent on the interactions that we 
computed. Above this energy scale a normal T^ like be- 
havior is recovered for the resistivity. Our results are in 
agreement with the observed resistivity in Bismuth. We 
examined several other consequences of the existence of 
such a mode and showed that it would lead to a double 
plasma edge in optical magnetotransport experiments. 
Recent measurements of the optical conductivity agree 
well with our predictions. 

The main contribution of our paper is thus to show 
the importance of such a plasmon mode in these sys- 
tems. This opens the way to many lines of work centered 
around the existence and role of such collective modes. 
First it is important to explore the consequences of these 
collective modes for other transport properties. One of 
the most important is of course the Nernst effect which in 
Bismuth is one of the largest reported. Analysis on how 
the plasmon modes can modify the Nernst transport is a 
non trivial questions that we plan to analyze in details. 
Another important direction is to go beyond the RPA 
analysis. This is necessary to investigate other semimet- 
als such as Sn-doped Bismuth or IT — TiSe2 or to take 
into account the effects of pressure on Bismuth. In that 
cases, the interactions between the plasmons play a more 
important role. We have suggested a phenomenological 
model which can potentially be used to tackle these ef- 
fects. However its analysis is highly non trivial and this 
will provide certainly an exciting line of investigations for 
the future. 



Appendix A: Acoustic plasmons in bismuth 

We give here some arguments and energy scales for the 
acoustic plasmons in Bismuth. 

The existence of acoustic plasmons was initially in- 
vestigated for Bismuth^* by means of a simple spher- 
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ical Fermi surface model with = ruh = 0.14 and 
nir — — 0.047. Based on this approximation it 
was concluded that the acoustic plasmon was very weak 
(fcc = 0.02kps), although the authors also mentioned in 
their conclusion the issue of Fermi surface elipticity. 

Solving when taking into account the full elipticity is 
a formidable, albeit necessary task, as is obvious when 
looking at the parameters given in Table HI The disper- 
sion relation of the collective mode is given by ^ in 
which the full masses should be inserted in the Il{q,uj) 
functions. Since the dispersion is quadratic one can in 
principle rescale the integration over fc in ([7]). One thus 
has 





R 


rv 


Cac 


-1 

T 


k* 


eRPA 


11 


0.2 


1.8 


0.4 


0.25 


UD 


17 


0.14 


2.1 


0.45 


0.3 



TABLE II. The parameters R and r„ (see text for defini- 
tion) of acoustic plasmon along the trigonal axis in the limit 
q — >■ 0. eRPA is an average taking partially into account the 
ellipsoidal shape of the Fermi surface, while UD is a purely 
unidimensional approximation. The good agreement between 
these two extreme cases suggests that the corresponding re- 
sults are probably representative of the experimental values 
in Bismuth 



(Al) 

where IP is the normal Lindhard function with all the 
masses set to 1. For the case oj/q — c constant and g — >■ 
the Lindhard function becomes particularly simple 



n"(g^ 0,a; = cq) 



kp 



-1 



log 



1 



1 



2n ^ 



(A2) 

where kp/2 — Ep the Fermi energy of the corresponding 
species. Note that the two Fermi energies are not neces- 
sarily the same for the two species of particles and must 
be determined by the neutrality condition. When going 
to higher order diagrams (vertex corrections) one real- 
izes that not all scattering directions are equally prob- 
able. This means that the coefficient in front of 11 in 
Ea lAll will change, because it depends on the anisotropic 
mass tensor. However the form of the Lindhard func- 
tion should not be affected by this procedure, thus the 
functional dependence of plasmon's dispersion relations 
should remain unchanged. 

One can thus try to simplify such an equation. An 
extension of RPA (cRPA) taking into account the ellip- 
soidal shape^l made the approximation to compute Cac of 
averaging the faster component on a spherical Fermi sur- 
face, while retaining the ellipsoidal shape of the heavier 
(slower) component. If one follows this procedure, and 
for example on the trigonal axis the electrons are faster 
than holes, so we need to substitute = m?i3 — 0.63 
which is already significantly larger. In the direction per- 
pendicular to trigonal the situation is even more complex. 
Holes are quite heavy along the trigonal axis direction 
which gives large spherical mass and suggests that they 
are always the slower component. However they are very 
light in the perpendicular direction (see Tab.|l]). It is then 
easy to forget that electrons, while extremely light along 
trigonal axis, are (even in average) significantly heavier 
than holes in any perpendicular direction. A spherical 
approximation (even the eRPA) would thus need to as- 
sumes to be able to do the average an opposite mass 
order than it is in reality! This suggest that a different 
approach is necessary. 



A simple average over a spherical Fermi surface comes 
from the summation over k present in I{^{q,uj) and 
the assumption that all scattering directions are equally 
probable. This is directly connected with neglecting in- 
teractions between electrons and holes within a bubble, 
and in particular the vertex correction within the bubble 
coming from the Coulomb interaction. Given the long 
range nature of the interaction, it is natural to expect 
that such corrections would enhance the small q scatter- 
ing. We can thus expect that the averaged mass that 
enters our expression for the dispersion relation is closer 
to the mass along the direction of q rather than the av- 
eraged mass. 

Although a full calculation is difficult and clearly be- 
yond the scope of the present paper we can have an idea 
of the importance of such effects by comparing to limit- 
ing situations, for the case of q along the trigonal axis. 
We take on one hand the averaging of the eRPA^ and on 
the other hand the unidimensional approximation where 
we simply retain the mass along q. For an isotropic 
case the expressions ([9]) considerably simplify and al- 
low to define two important parameters for the equa- 
tion namely the relative weight between the two terms 
R = {nx_p I yi^-pY = i'm-skps) / (nirkpr) (where Klp_p 
was defined in the context of Eq l20|) and r„ = Vps/Vpr 
which is the ratio of the two velocities. The results are 
given in Table. [Hi We see that two extreme approxima- 
tions give rather similar values. The correct values for 
Bismuth should therefore be between these two extremes. 
This suggests that k* is far from being negligible. Tak- 
ing such values for fc* gives an energy of T* ss QK for the 
temperature below which the acoustic plasmon becomes 
strongly coherent. 

In the perpendicular plane the situation is significantly 
more complex. In fact we have to work with a multi- 
plasma problem given the three electron pockets so we 
can only give qualitative arguments. Let us stay within 
the UD approximation, remembering that it tends to 
slightly overestimate k* . From Fig. [1] we immediately 
see that by changing the q direction we are facing the 
following situation: 

• there is one central hole pocket with an angle inde- 

1 2 

pendent mass m, ' = 0.067. 
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• the mass of each electron pocket can vary drasti- 
cally, from very light — 0.0015 up to the heav- 
iest of all ml = 0.198 

• the electron pockets are situated with an relative 
orientation of 7r/3 and thus for every q direction 
there are either very light and/or very heavy elec- 
trons present. 

This open the possibility for several plasmons to emerge, 
but most of them are strongly Landau damped. A rough 
estimation gives (for q not far from the bisectrix between 
two pockets): for light hole-heavy electron i? « 1.15 
k* w 0.07; for light electron-heavy hole i? « 1.2 ^ fc* w 
0.1; for light electron- heavy electron i? w 13 ^ fc* w 
0.22. For another intermediate orientation (tt/G angle 
with the bisectrix) both electron masses will be smaller. 
A more precise evaluation of k* would be quite difhcult 
and probably also susceptible to a significant error. Let 
us only emphasize that for the Fermi energies Ep sa "^Ep, 
thus a smaller value of k* in the perpendicular plane does 
not immediately imply a smaller T* in this plane. 

Finally let us make two comments for plasmons far 
from the trigonal axis: 

1. If we had taken the cRFA approximation we would 
have found rather low values of fc* along the bi- 
sectrix axis which would imply a notable trigo- 
nal/bisectrix anisotropy of the observed T*. This 
is not the case found experimentally. 

2. Due to particular positioning of Fermi surface pock- 
ets it is practically always possible to find carriers 
with rather different masses. Thus we expect that 
a drop of fc* for some directions is rather the ex- 
ception than the rule. These rare cases should not 
affect significantly the q momentum integrals such 
as the one in (|4T]). 



Note that in the above estimation of the T* values 
we have neglected local field corrections. These always 
tend to reduce T***^. Precisely: the larger fc* the larger 
influence of local field corrections. The "Deybe" temper- 
atures found by resistivity fits by UherJA are thus lower 
but this possible discrepancy was already discussed in the 
Secinill 



Appendix B: Anisotropy of the Baber resistivity 

In order to compute the anisotropy of the resistivity 
due to the Baber scattering we start from the formula 
obtained in Ref. '27' (Eq. 2.14 in this paper) and apply it 
to our problem with slow and rapid carriers: 

= 4kJiUkm - fiUk) + A,)) 

fiUpm'fi&ip)-^)) (Bi) 



where Ag is the energy exchanged in the scattering pro- 
cess (determined by the exchanged momenta q) and 
Pq = TTj'^/{nDx), Dx being the density of carriers in 
a given direction, J the strength of carrier's interac- 
tion. In our model, instead of the J of (jBip . included in 
Pq, the averaged scattering rate W must be substituted. 
Following the derivation of Ref. [23 we take the carri- 
ers velocity along the resistivity axis as the variational 
parameter = Vx- According to textbook procedure 
(p. 283 in Refl43l) Vpkq is then equal to the total change 
of $fci during scattering event which in our case means 
Vpkq — AVx- When one derives the relation AVx{q), one 
needs to take the masses along x-axis. The momentum 
along a given direction must be conserved, which implies 
that |AT4((7)| ^ QxlmJx ~ "^hi I where q^ is the momen- 
tum exchanged along resistivity axis. 

The summation in (jBl|) is taken over vectors of incom- 
ing and outgoing carriers. The detailed, all-T treatment 
of (jBl|) would be quite complicated, but we are interested 
only in lowest temperatures, where the energy exchanged 
during scattering goes to zero A — >■ (we always stay on 
the Fermi surface) and the density of both electrons and 
holes is constant (the chemical potential in semimetal is 
rather stable). We also know from band structure that 
there is no nesting fcg 7^ k^ and these Fermi wavevectors 
are also constant. 

Note that the derivation in Ref. [13 was done for a 
two dimensional system. Working in three dimensions 
requires one extra angle to complete the spherical co- 
ordinates. It is an angle between incoming momenta of 
light and heavy carriers between Z(fc,p) which was ob- 
viously equal to zero for 2D case. As long as we are 
interested in the resistivity along the high symmetry di- 
rections (which is the one usually measured) and in the 
limit of small A we may assume that the functional re- 
lations A{dj:^,9p^) and \q\{dj:^,dpq'), originating from en- 
ergy conservation, do not depend on the new angle 
One integrates out this extra variable, but this does not 
change the low energy T-dependence (because the con- 
dition A ^ allows to give a unique relation between 

Oj:^,dpq)- 

With the above remarks we can evaluate the mo- 
mentum sums, transposed into integral over the angles 
^fc(f' %9' ^- masses rus (0^^) and rrir {Opq) entering 

to formulas accounting for energy conservation (during 
scattering event) A{9^^,6ff^) and kK^'^:^, %q) are aver- 
aged over all angles on the Fermi surface when integra- 
tion over all possible orientations in taken. This brings 
us to an analog of Eq. 3.4 in Ref. '27 with the anisotropy 
factor extracted: 



Px 



(AVx/qx 



J 



(cos(6l£g-) - cos(6ip-,-))^ 



4sinh^(A,-;3/2) 



(B2) 



where an integral is angle averaged which also means 
averaged over all momenta on the Fermi surface. The 



17 



integral is rather complicated in 3D. We can safely as- 
sume that the anisotropy of the resistivity is caused 
only by the two factors in front: a change of carri- 
ers velocity AVx (along the resistivity direction) during 
the scattering event which we already discussed above 
and the denominator Dx- The denominator Dx is con- 
stant when fJ.{T) = cste and it introduces the proba- 
bility of such a scattering event, which is proportional 
to carrier density determined along the given direction 
by Dx ~ {n^ex + ™/^a;)^- The resistivity can be ex- 
pressed as px ^ \AVx\/Dx. This implies that, in the 
low temperatures (when p ^ T^), the main contribution 
to the anisotropy of resistivity is accounted by the fac- 
tor \AVx\/Dx which after straightforward simplification 
gives: 



1 2 



(me 



rrihi 



(B3) 



where / is angle averaged integral. In the case of two 
pockets Fermi surface this immediately leads to: 



Px 
Pz 



ruhxl/imex + ruhx) 



(B4) 



'riez - mhz\/{mez + ruhz) . 
rUe and m/j are second rank tensor with values (along 



main axis) given in Tab. HI A reminder is necessary: 
the masses which are substituted into (f2T|) are not the 
ones given in Tab. ID Indeed the masses in (|2ip are the 
ones along external directions while the Tab. U gives the 
ellipsoid parameters. This point is particularly important 
for electrons along the trigonal axis: there is very small 
mse but (due to out of bisectrix plane tilt •§ of the electron 
ellipsoid) vriez contains an admixture of relatively large 



mie-. 



-/=m3-/cos(^)2+mr/sin(^)^ 



(B5) 



There is one extra complication, not present only if the 
selected direction is the trigonal one. It comes from the 
fact that in Bismuth we have more than two families 
of carriers. In the particular case of the bisectrix axis 
(along which resistivity is frequently measured) we have 
to deal with three different types of carriers (holes, light 
and heavy electrons) with certain probabilities for each 
type of scattering event. This leads to: 

3 , , 6 , ^ 2 , 

Px = YlPieh(mei,mhi) +—p2eh[me2,rnhi)+—pee{rnei,m 

(B6) 

which one can compare with pz measured along trigonal 
axis precisely what we did in Sec. Ill CI 
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